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Abstract 

We propose a new approach towards the understanding of confinement. Starting from 
an anisotropic five dimensional pure gauge theory, we approach a second order phase 
transition where the system reduces dimensionally. Dimensional reduction is realized via 
localization of the gauge and scalar degrees of freedom on four dimensional branes. The 
gauge coupling deriving from the brane Wilson loop observable runs like an asymptotically 
free coupling at short distance, while it exhibits clear signs of string formation at long 
distance. The regularization used is the lattice. We take the continuum limit by keeping 
the ratio of the lattice spacing in the brane over the lattice spacing along the extra 
dimension constant and smaller than one. 
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1 Introduction 



Four dimensional SU(N) gauge theories have a unique phase in which interactions are 
confined [1]. The only known fixed point in their phase diagram which is at vanishing 
gauge coupling, g 2 = 0, is an ultraviolet (UV) fixed point where weak coupling perturba- 
tion theory is a good description. According to perturbation theory the interactions in 
this regime are dominantly of Coulomb type, with the charge obeying the renormalization 
group (RG) flow appropriate to an asymptotically free coupling [2]. Perturbation theory 
is however oblivious to the long distance effects of confinement. In order to see those, 
one must increase the coupling to larger values where the only probe we know of are 
lattice [3] Monte Carlo (MC) simulations. 1 As the MC simulations reveal, globally the 
static potential seems to consist of two distinct regimes. At short distance it is indeed 
of a Coulomb form but at long distance it develops a linearly growing behavior. More 
precisely, at short distance the coupling a = g 2 /4n defined as (for SU(2)) 

««(l/r) = ~\c{r) , (1.1) 

with c(r) = (l/2)r 3 F'(r) and F the static force, decreases according to the perturbative 
RG flow and around a certain scale r s , defined as the scale where perturbation theory 
breaks down, c(r) plateaus around the value —0.3 [4]. The physically motivated expla- 
nation of this behavior comes from an effective string description according to which 
confinement results in the formation of a string like flux tube. The positive slope of 
the linear term in the static potential is then interpreted as the tension a of this string. 
The massless degrees of freedom that describe the fluctuation of the tube are Goldstone 
modes with an effective action that can be written in the form of a derivative expansion 
constrained by Poincare (and perhaps also diffeomorphism) invariance. The kinetic term 
in this class of "world-sheet" actions, when integrated out, yields a universal term, the 
Liischer term [5], with value — (d — 2)7r/(24r) in ci-dimensions. This is the term that the 
MC sees as a plateau in c(r) at large distance. Similar properties are believed to hold for 
any generalization of the pure gauge theory where the string is stable [6] . In summary, the 
description of the static potential from short to long distance entails two different analytic 
methods, namely weak coupling perturbation theory and an effective string description, 
each of which is blind to the physical effects that the other sees. The two approaches are 
bridged by the MC which, in principle, can probe the whole static potential from weak 
to strong coupling. In practice however [7], at short distance, simulations tend to give 
a much less precise description compared to the usual continuum field theory Feynman 

1 At very strong coupling we can also use strong coupling expansion methods but the information we 
can extract from them is limited and in addition there is no analytic connection to the weak coupling 
regime. 
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diagrams due to the enhanced lattice artifacts and the same applies to any analytic per- 
turbative lattice computation. It is fair to say that a global analytic understanding of the 
static potential of confined 4d gauge theories is missing. 

In [8] we proposed a regularization of four dimensional gauge theories which could allow 
for such a unified description. The idea is to start form an SU (N) pure gauge theory in 
5d and via a systematic expansion around an anisotropic mean-field background [9] try 
to reach an ultraviolet fixed point (or points) in the interior of the phase diagram where 
a second order phase transition takes place. The anisotropic background has the value 
TJo5 along the extra dimension and vq along the other four directions (see Appendix A). 
In the confined phase it vanishes while in the deconfined phase it is non-zero along all 
directions. There is a (unstable according to the leading order mean-field method) phase 
where Vq ^ and 1% = 0, the layered phase. Only on the isotropic lattice the background 
is isotropic. A line of second order phase transitions was found on the boundary between 
the deconfined phase and the layered phase. From the side of the deconfined phase, near 
the phase transition, surprisingly, the system reduces dimensionally. Even though not 
the same, the physics in this phase is similar in spirit to the layered phase of [10]. We 
have called in [8] this phase the "d-compact phase". The low energy degrees of freedom 
of the dimensionally reduced system are those of the four dimensional Georgi-Glashow 
model which is in the class of theories described in the first paragraph. The technical 
tool used to carry out the expansion is the lattice regularization. This allows for a well 
defined description and control of the quantum theory but there is a price. Because the 
final results for the observables are expressed in terms of finite lattice sums it is not easy 
to take the continuum limit analytically. We emphasize though that this is merely a 
technical obstruction. The method being fully analytic should, in principle, allow one to 
carry out all limits without resorting to numerical methods. This will be attempted at a 
later stage. Here we perform a numerical analysis of the analytic results of [8], with our 
main focus on the static potential. We will take carefully the continuum limit and we will 
try to argue that the static potential oriented along the four dimensional hyperplanes, 
computed in this scheme, reflects both the asymptotically free and the confining aspects 
of the coupling. 

The mean-field expansion comes with certain caveats. The mean-field background is 
gauge dependent. Within the class of Lorentz gauges we found our physical observables 
independent of the gauge fixing parameter £ to leading order [8]. There is no guarantee 
that the expansion converges. It is known however that the corrections come multiplied 
by powers otl/d and therefore convergence is expected to become better as the number of 
dimensions d increases. The mean-field sometimes fakes phase transitions. Even though 
the known such cases are generally less sophisticated compared to our construction, it is 
still conceivable that there is an intricate way in which the mean-field does generate a fake 
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picture. This is where the MC investigation of the phase diagram of these theories will be 
crucial [11]. The reason we proceed with presenting the results of the mean-field ignoring 
this possibility is that even though the MC may not see an UV fixed point where the 
system reduces dimensionally, the approach seems to have an independent value. It can 
serve as a new analytic laboratory of confining gauge theories from short to long distance. 

2 The model and its Lines of Constant Physics 

The model considered in [8] is an anisotropic SU(2) lattice gauge theory in 5 dimensions 
defined in a mean-field background. Gauge theories in five dimensions are defined on 
an anisotropic, infinite, hypercubic lattice via two independent dimensionless parameters. 
One is (3, the lattice coupling and the other is 7, the anisotropy parameter. One way 
to define them is through the dimensionful parameters of the lattice. We consider first 
finite hypercubic lattices with the same number of lattice points along the four dimensions 
L = l/di and along the fifth dimension L = 2nR/a 5 and eventually take the L — >■ 00 limit. 
/ and R are the physical sizes and as appropriate to anisotropic lattices we take different 
lattice spacings along the four and fifth dimensions. The coupling of a five dimensional 
gauge theory is denoted by g§ and it has dimension of ylength. The anisotropy parameter 
can be defined at Oth order as 7 = a 4 /a 5 = l/(2nR) and the SU (2) coupling as (3 = Aa^/g^. 
In these variables, the perturbative regime is located at /3 — 00. The above mentioned d- 
compact phase appears instead for (3 ~ 0(1) and 7 < 1, obviously far from perturbation 
theory. The line of second order phase transition extends in a range that corresponds 
to approximately 7 < 0.62. Physically, this is a situation where the extra dimension is 
larger than the spatial directions and gauge interactions are localized on four dimensional 
hyperplanes. 

The action used to compute observables was the Wilson plaquette action 



where the first term contains the effect of all plaquettes U p along the four dimensional 
slices of the five dimensional space and the second term contains the effect of plaquettes 
having two of their sides along the extra dimension. In the mean- field, for SU(2), the 
fluctuating degrees of freedom are complex valued quantities V(n, M) located at the lattice 
site n and pointing in the direction M. The schematic expression for the correction to 
the expectation value of a physical, gauge invariant observable O to second order in the 




(2.1) 
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mean-field expansion has the form 



<o> = 



+ 




i,k,l,m 



(2.2) 



where K is the appropriate lattice propagator and the sums over i,k,l,m are sums over 
links. Derivatives and contractions are taken in the mean-field background. All such 
quantities are gauge independent. Polyakov loops with scalar and vector transformation 
properties represent corresponding classes of states in the Hilbert space and from the 
exponential decay of their Euclidean time correlators their mass spectra can be extracted. 
The expectation value of Wilson loops can be used to extract the static potential. The 
anisotropic lattice admits two inequivalent classes of Wilson loops. One class consists 
of the loops oriented along the time direction and one of the spatial directions. The 
other class consists of loops along the time and the extra dimension, the latter defined 
as the direction along which the background is different from the background along the 
other four directions. We call the static potentials derived from these two classes as V4 
and V5 respectively and define as Fj(r — %/2) = {Vj(r) — Vj(r — Q>j)}/o>j, j = 4,5 the 
corresponding forces. 

The dimensionless vector mass (see Eq. (A. 15)) is found to have a dependence only 
on the lattice size [8] 



where cl is a constant with numerical value 12.61. Therefore the system can not be in 
the Higgs phase in a mean-field background. Hence, if dimensionally reduced, it must 
be in a confined phase. The scalar mass a 4 m# (see Eq. (A. 12)) can be used to measure 
a critical exponent of v — 1/2 as the second order phase transition is approached. The 
potential Vs(r) (see Eq. (A. 10)) determines the interaction of static quarks along the extra 
dimension. The static potential V^(r) (see Eq. (A. 9)) is the quantity that dictates the 
behavior of the dimensionally reduced coupling and it is well defined for any value of the 
distance r. For more details on these computations we refer the reader to [8]. 

It is convenient to trade the bare lattice coupling [5 for the dimensionless physical ratio 
p = a^mw I aimH and parametrize the model by the two dimensionless physical ratios p 
and 7. The lattice spacing can be then chosen to be measured by a^mn- An issue arises 
when one realizes that on an infinite lattice both the numerator and the denominator 
of p approach zero in the limit 04 — > 0. This means that the continuum limit must be 
reached in a way that regularizes p. The method to achieve this is to approach the phase 



a±m w = c L /L , 



(2.3) 
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Figure 1: Left Plot: the function H{x) in Eq. (3.3) along the LCP A. Right plot: the ratio of 
the extra dimensional to four dimensional force F^{r) / F^{r) along the LCP A. 



transition at fixed 7 and at every step, adjust L and (3 so that p remains constant. When 
taken in this way, the infinite lattice and continuum limits coincide: 

(L — > 00, f3 — > /3 c )l7,p=const •<=>■ continuum limit (2.4) 

A crucial fact that allows one to apply this method without obstacles is the fact that 
a^rriw depends only on L and that a 4 m# depends only on (3 and 7. We call from now on 
such continuum limit trajectories on the phase diagram approaching the phase transition 
Lines of Constant Physics (LCPs). The whole discussion above can be generalized to any 
SU(N) gauge group. Here we concentrate only on 577(2). 

3 The continuum limit 

Following [12] we define a physical scale through the condition 

r s 2 F 4 (r s ) = s = 0.2. (3.1) 

s = 0.2 is chosen so that r s ja^ lies in the transition region from the short distance behavior 
of the force to the long distance one. We will first consider the LCP trajectory defined as 

LCP A : 7 = 0.55 and p = 0.625 , (3.2) 

then consider other values of p for the same 7 and finally we will analyze other LCPs 
labelled by values in the range 0.62 > 7 > (i.e. move down along the line of second 
order phase transitions). 

In the left plot of Fig. 1 we show the dimensionless function 

H{x) = F 4 (ry\ r=xrs (3.3) 
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for various lattice spacings corresponding to L values in the range 100, . . . , 400. The 
function shows good scaling. The shape of H(x) is different than in a pure SU(2) gauge 
theory, see [12], in particular it decreases as x increases. We will see below that this is 
due to the presence of a positive 1/r 2 contribution at short disctance and a large negative 
logarithmic contribution at large distance in the static potential V4. 

In the right plot of Fig. 1 we show the ratio F 5 (r)/F 4 (r) on a logarithmic scale as 
a function of r/r s for different lattices along the LCP A. When we take the continuum 
limit the ratio decreases steadily which implies that the force F 5 is not physical. This 
is an evidence for localization: since F 4 has a finite continuum limit (see left plot in 
Fig. 1), F 5 tends to zero, implying that the five dimensional space decomposes into a 
set of non-interacting Euclidean four-branes along the fifth dimension. On each of these 
branes the localized light degrees of freedom are therefore expected to be those of the four 
dimensional Georgi-Glashow model. Since the mean-field background does not break the 
gauge symmetry spontaneously, these degrees of freedom must be in the confined phase. 
In the following, we will take the continuum limit of the static potential V4 and try to 
argue that our picture of dimensional reduction is consistent. 

To extract physical information from the static potential Eq. (A. 9) we start from an 
ansatz of the form 

v 

with p > — 1 (p = corresponds to a logarithmic term and c_i = a), solve locally for the 
coefficients c p and plot them as a function of r (in practice we will rescale all dimensionful 
quantities by appropriate powers of r s to make them dimensionless) . The limit L — > 00 
will be taken by computing the coefficients for increasing values of L and extrapolating 
to the infinite lattice. 



3.1 Short distance 



For r/r s < 1 we choose the ansatz 

V(t) = vl + ^ + ?1, (3.5) 

the role of the 1/r 2 term being the check of an imperfect dimensional reduction. In Fig. 2 
we compute c\ and ci (see Appendix B) and compare the coefficient c\ to the 1-loop and 
2-loop RG evolution formula 

A w -r=(M 2 )- 6l/(26i) e- 1/(2b092) , (3.6) 

where g 2 = — (167r/3) C\ and g 2 = g 2 (l/r). For the Georgi-Glashow model b = 7/(Att) 2 
and 61 = (115/3)/(4tt) 4 . The 1-loop formula is obtained for 61 = 0. Unlike for QCD 
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Figure 2: The short distance coefficients c p for the LCP A. 



there is no available experimental or lattice MC prediction for a scale, so we will choose 
it according to convenience to A q qr s = 0.277. Also, due to the high degree of the discrete 
derivatives involved in the determination of Ci(r), the latter can be defined only for 
r/<24 > 2 which makes the comparison with the RG formula valid near r = hard. 
Nevertheless, it is clear from the left of Fig. 2 that even though the continuum limit has 
not been reached, as L increases there is a stable convergence to the theoretical curve. 
The continuous line on the left is the 1-loop curve and the dashed one on the right is 
the 2-loop curve. To illustrate our point we have shifted the 2-loop curve by a constant, 
shown as a dashed dotted line. In addition, according to the right of Fig. 2, as L — > oo, 
the 1/r 2 piece gradually disappears. 



3.2 Long distance 

For r/r s > 1 we choose the ansatz 

V(r) = /i + ar + c log(r) + - + ^| (3.7) 

and compute the coefficients using the discretized formulae in Appendix B. The plots 
show plateaus forming for all four coefficients at the same range of distances 

r/r s e [2.15, 2.80], (3.8) 

a sign that the ansatz is close to the actual form. Moreover we have checked that the 
plateaus form essentially independently from the way we discretize the derivatives of V(r). 
We take as the plateau value of a coefficient the average of the quantity over the range 
indicated in Eq. (3.8). Reading off the value of a coefficient (or of a mass) from a plateau 
introduces errors, which we add on its plot. In Fig. 3 we compute the continuum limit of 
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Figure 3: The continuum limits of the long distance coefficients c p for the LCP A. 

the c p for the LCP A using a linear fit in a\, which is the expected form of leading cut-off 
effects in the static potential [7]. The continuum limit value of the string tension is 



ar s = 0.0433(13) 



(3.9) 



and using the scale K q q estimated in Section 3.1 we obtain roughly A qq -/y/a ~ 6.4. Even 
though not from the same theory, the pure four dimensional SU(2) gauge theory value 



[12,13] (which is the only quantitatively reliable case we know) A^ M) /v^ - 0.68(10) tells 
us that the string tension is probably rather small in our model. A physical understanding 
of this fact can be obtained by reading off the continuum limit value 



y/al = 0.69 



(3.10) 



from Fig. 3 and Fig. 4: the box seems to be barely sufficiently large to sustain a stable 
string. In order to increase this number, one would need to increase I, which as we will 
argue in the following, requires decreasing p and increasing L. On the other hand, the 
continuum limit value of the coefficient C\ is 



ci = -0.2586(23) , 



(3.11) 
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Figure 4: The continuum limits of r s rri}j and l/r s . The errors of nifj originate from the ones of 
mw and come from the uncertainty in the effective mass at L = 100. mjy is scaled proportional 
to 1/L. On the right linear fits in a\ show the continuum limit. 



in excellent agreement with the universal value of the d = 4 Liischer coefficient —n/12, as 
it should. The somewhat surprising term that seems to be though necessary for a correct 
interpretation of our formulae is the log(r) term. In fact, if left out, a consistent picture 
can be hardly obtained. Regarding the 1/r 2 term that we seem to be also seeing, as well 
as other possible higher negative powers of r we would not like to make any committing 
statements until a more complete understanding of the structure of the lattice formula 
Eq. (A. 9) both from the field theory and the effective string point of view is obtained. 

The next issue concerns the universality of the continuum limit. It is easier to formu- 
late the question in physical terms. In Fig. 4 we compute the continuum limits of r s rriH 
(left) and of l/r s (right) for three RG trajectories at 7 = 0.55 with different p values of 
0.4, 0.625 and 1. Evidently, r s m# has (within errors) a p independent continuum limt, 
around 1.25. p on the other hand, when L is varied, controls the physical size of the 
lattice, i.e. the size of the box: as p increases, the box shrinks. When L is kept constant, 
it controls the lattice spacing: as p increases, the lattice spacing decreases. It is interest- 
ing to notice that since p {l/r s ) has a universal continuum limit which is around 10, and 
because (by definition) p{l/r s ) {r s m H ) = c^, the continuum limit value of cl is, within 
its error, the same as its finite L value 12.61. Finally the critical law of approach to the 
phase transition along an LCP can be expressed as [8] 

18 = A 7 -2, (3.12) 

with k some constant, confirming that it is indeed the quantity pL which controls the 
limit. 

To approach the continuum limit at a given value of 7 (at a sufficiently large L) it is 
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convenient to first tune p to an optimal value and then take L as large as possible, the 
latter being the more time consuming step. The larger 7 is, that is the more isotropic the 
box becomes, the smaller box is appropriate to sustain a string (in a large box a small 
anisotropy is nearly unnoticeable), which means that a larger p will be needed. But the 
box can not be arbitrarily small. The string being a long distance effect, will be unstable 
in a too small box. Thus, we expect to have an upper limit on p 

P < Pinst (at given L) , (3.13) 

beyond which some sign of instability should appear. Indeed, beyond a certain p the 
string tension turns negative, a typical sign of instability. On the other hand when 7 is 
small, the physical volume of the four dimensional branes shrinks compared to the extra 
dimension and therefore a large box, i.e. a small p is necessary to sustain a string. But 
the box can not be arbitrarily large either. A lower bound 

p > p dr (at given L) (3-14) 

on the possible values of p exists, beyond which the critical law Eq. (3.12) is not valid 
anymore and dimensional reduction is lost. The coefficients a, cq, c\ plateau at zero with 
only C2 remaining non-trivial, as appropriate to a five dimensional static potential. Also 
the ratio F 5 /F 4 approaches a constant, which is 1 for the isotropic case. 

The above discussion implies essentially that only a box whose physical size and 
anisotropy are within a certain allowed (correlated) range can describe a continuum four 
dimensional theory with a stable string. It turns out that in our case this range for the 
anisotropy parameter is around 

0.6>7>0.5, (3.15) 

a rather narrow window. Our prime example, LCP A (for which 7 = 0.55), falls in the 
middle of this range. As either of the two bounds in Eq. (3.15) is approached the system 
seems to require very quickly huge lattices in order that the continuum limit is reached. 
For 7 = 0.55 the physical size of the lattice can be approximately in the range 2 > p > 0.5. 
Outside this range we observe either a negative string tension or five dimensional physics. 
Again, soon deviating from the LCP A value (p = 0.625), the continuum limit seems to 
demand larger lattices than we can afford at the moment. 

The remaining question is how to construct such a box that fits the "universe" . The 
answer is in Fig. 4 and Eq. (3.12): we start with values p\ and L\ that describe four 
dimensional physics (for LCP A this is realized at p x = 0.625 and L x = 400) and decrease 
p 2 < pi adjusting the lattice size so that L 2 > (pi/ p 2 )L 1 . Like this f3 in Eq. (3.12) does 
not increase and the bound Eq. (3.14) is not violated. In particular these operations 
should make it possible to eventually reach y/al > 1. Finally, it is interesting to observe 
that in this limit the scalar mass becomes much larger than the vector mass. 
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4 Conclusions 

A five dimensional pure gauge theory on an anisotropic lattice can be described by fluc- 
tuations around a mean-field background. The phase diagram has a line of ultraviolet 
fixed points where the system can reduce to a collection of non-interacting Euclidean 
four-branes and the continuum limit can be taken. The static potential together with the 
masses of the lightest fields can be used to describe the system away from perturbation 
theory. The four dimensional gauge coupling derived from the static potential at short 
distance runs like an asymptotically free coupling while at long distance it becomes the 
Liischer coefficient of a confining string. The four dimensional theory recovered in the 
continuum limit is not a pure SU (2) gauge theory. The potential shows confinement but 
also has a large logarithmic contribution and the vector particle is lighter than the scalar. 
This theory should correspond to a region in the parameter space of the Georgi-Glashow 
model. 

Even though we have computed all the observables analytically, we did the continuum 
limit analysis numerically. Clearly, the success with which the model passed several severe 
tests calls for further study, where the lattice propagator is inverted analytically and the 
finite lattice sums are performed explicitly. Then the continuum limit itself could be 
studied analytically. 

The static potential is computed using a MATLAB code. Using an AMD Phenom 
II 3.2 GHz processor, the computational cost is 35 core-hours on a L = 200 lattice and 
about 34 core-days on a L = 400 lattice. 

Acknowledgments. We thank R. Sommer for correspondence. N. I. acknowledges 
the support of the Alexander von Humboldt Foundation via a Fellowship for Experienced 
Researcher. 



A Appendix 

In this Appendix we summarize the expressions for the observables computed analytically 
in [8] and analyzed numerically here in the main text. 

The propagator in momentum space is defined as 
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(A.4) 

The above propagator is written in the Lorentz gauge with parameter £. Also, we have 
used the following notations and conventions: The propagator is an object with the 
following index structure: 

K(p',M',a';p",M",a") (A.5) 



with p the discrete momentum 



k M = 0, 1, . . . ,L - 1, 



(A.6) 



M = 0,1, 2, 3, 5a Euclidean index and a = 0, 1, 2, 3 an index taking value in the Lie group 
SU(2). The Euclidean structure is built in the matrix form. We use the notation sm = 
sin (p' M ), c M = cos (p M ), s M/2 = sin (p' M /2), c M/2 = cos (p' M /2). The only special case not 
explicitly shown in the matrices is that in the diagonal elements, c 5 = 7 2 ^r cos (p' 5 ). On 
the diagonals Y' implies summation over all Euclidean indices leaving out the one that 
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corresponds to the row/column index of the term. The functions h,2,3 are the usual Bessel 
functions of the type /. The background values v$ and U05 along the four dimensional 
branes and extra dimension respectively, are determined by the extrimization of 

" H 1 ' H ! j 'll^ l-Ul-\)T;ri :i ■ (d- 1 ) / /(77 ) • „(7^) • (d l)h v +h 05 Vo 5 



L 5 



7 



(A.7) 



which yields the conditions (primes denote derivatives with respect to the argument) 



— — 6/3 

v = -u(h )' , ho = — v 3 + 2/3-fv vl 5 , 

7 

^05 = -u(h 05 y , h 5 = 8/?7^05 , 



(A.8) 



for d — 5. In the above we have introduced the function u(x) = — In . On an 

anisotropic lattice there are two inequivalent Wilson loops of spatial length r. One along 
the four dimensional hyperplanes for which the static potential is given by 



V 4 (r) = -21og(u ) - x 



2v 2 L* 



E 



Pm^O'Po— 



^(2cos(p' fe r) + 2) 



+3 E 



]T(2cos(p' fc r)-2) 



CooO/,!) 



(A.9) 



The one along the extra dimension is given by 
V 5 (r) = -2\og(v ) 

~ ^77 E { [2 cos(p' 5 r) + 2] C m \p\ 0) + 3 [2 cos(^r) - 2] ^-^ ^ 



(A.10) 



The matrix C in Eq. (A.9) and Eq. (A. 10) is defined from 

K(p',M',a';p",M",a") = <W'<WCW»(p', «')■ 
The scalar mass is derived from the correlator 

Ch(*) = ^rf^E^^oOEl^^)! 2 ^ 1 ^^,^^^; (pUj/ 6 ),5,0) > 



(A.H) 



(A.12) 

where Pq ^ is the associated Polyakov loop evaluated on the background (being an overall 
constant, its value is irrelevant for the mass) and 



L ~ 1 ip(m+l/2) 



A(L) w = E^ 



m=0 
14 



^05 



(A.13) 



For the W gauge boson mass we define the contractions 
and 

T~\t,f?,a) = J2^' ot K~\(PoJ),5,a). (A.14) 

P'o 

in terms of which the vector correlator reads 

CVW = ^(^o (0) ) 4 N(0)) 4 ^^sin 2 (p' fe ) (T\t,i?,l)) 2 . (A.15) 

tf k 

The lightest state's masses are read off to second order from the exponential decay of the 
correlators according to 

(7(1) (t) + C {2) (t) 
m "^ ln ^)(t-l)+^)(t-l) - (A - 16) 
The scalar mass is non-trivial at first order while the vector mass is non-trivial at second 
order. 



B Appendix 

Here we write the discretization formulae for the computation of the local coefficients of 
Eq. (3.5) and Eq. (3.7) to fit the potential along the four dimensional hyperplanes. We in- 
troduce the dimensionless potential a^V(r) = V(x), define x = r/a^ and the dimensionless 



coefficients a = a\a, c 2 = c 2 /a 4 and c = a 4 c . C\ is already dimensionless. 
We use derivatives of the potential improved to 0(a|): 

V\x) = {V(x + l)-V(x-l)}/2, (B.l) 

V"{x) = V(x + l) + V(x-l)-2V(x), (B.2) 

V"\x) = (V(x + 2) -V{x-2) -2[V(x + l) -V(x- l)]}/2, (B.3) 

V""{x) = {V(x + 2) + V(x-2)-A[V(x + l) + V(x-l)} + QV(x)}. (B.4) 

In terms of these derivatives the coefficients in Eq. (3.5) are estimated to be 

Cl (x) = 2x 3 V'\x) + l/2x 4 v"\x), 

~ 2 {x) = -l/2x A v'\x) - l/6x 5 V"\x) , (B.5) 

and those of Eq. (3.7) 

a(x) = V\x) + 3xV"(x) + 3/2x 2 V"\x) + l/Qx 3 V""(x) , (B.6) 

ct(x) = -6x 2 v"(x)-4x 3 v'"(x)-l/2xW'"\x), (B.7) 

Cl (x) = -Ax 3 v"(x)-7/2x 4 v"\x)-l/2x 5 v"'\x), (B.8) 

ci{x) = l/2x 4 v'\x) + l/2x 5 v"\x) + l/12x 6 V'"\x). (B.9) 
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